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Abstract 

The internal electric field of an illuminated liquid droplet 
is studied in detail using both wave theory and ray theory. The 
internal field attains its maximum values on the caustics within 
the droplet. Ray theory is used to determine the equations of 
these caustics and the density of rays on them. The Debye series 
expansion of the interior field Mie amplitudes is used to calculate 
the wave theory version of these caustics. The physical 
interpretation of the sources of stimulated Raman scattering and 
fluorescence emission within a liquid droplet is then given. 

Introduction 

When a spherical liquid droplet is illuminated by a light 
source, the light is partially focused inside the droplet, forming 
patterns of caustic lines. The purpose of this paper is to analyze 
these patterns using a reformulation of Mie theory. One limitation 
of Mie theory is that it often gives little physical insight as to 
the origins of various scattering effects. As an alternative, ray 
theory often provides physical insight into the scattering process. 
However, ray calculations are incomplete and do not take into 
account wavelike phenomena such as diffraction and interference. 
The gap between Mie theory and ray theory can be bridged with the 
use of the Debye series expansion. The Debye series rearranges the 
Mie equations into an infinite series of internal reflection and 
transmission terms. Since the Debye series is mathematically 
equivalent to the Mie equations, it is exact and accounts for all 
the mechanisms that produce scattering. Its usefulness comes from 
calculating the individual terms of the series because they are 
analogous to reflected and transmitted light rays. These terms of 
the Debye series provide the physical insight into the origin of 
the caustic patterns inside the spherical droplet. 

This paper represents the first time to our knowledge that the 
Debye series has been used to study in detail the interior 
scattering pattern of illuminated liquid droplets. This paper 
also makes use of the full Mie theory, as well as ray theory to 
analyze the internal caustic structure of a liquid droplet. The 


results of stimulated Raman scattering (SRS) and fluorescence 
emission experiments are investigated as well, and the physical 
mechanisms that produce the sources of SRS and fluorescence in the 
droplet are determined using the above analytical methods. 

Ray Theory of the Caustics within a Liquid Droplet 

In this section we consider a spherical droplet of radius a 
and refractive index n whose center is at the origin of the 
coordinate system and which is illuminated by a family of rays 
propagating parallel to the positive Z axis. This geometry is 
shown in Fig. 1. Consider a ray whose angle of incidence at the 
droplet is 6,. The fraction of this ray's intensity that is 
transmitted into the droplet is T(0;) and the fraction that is 
reflected at the surface is R(0 ; ). For a>>\, the droplet surface 
is approximately flat, thus T(0;) and R(6,) are given by the Fresnel 
coefficients. 1 For the ray that is transmitted into the droplet, 
the fraction of its intensity that is internally reflected at each 
subsequent interaction with the droplet surface is also R(0j). 

We define the p-ray family to be those rays within the droplet 
for all that have made p - 1 internal reflections. For example, 
when the incident rays are transmitted into the droplet, they 
become the p = 1 rays, after their first reflection inside the 
droplet they become the p = 2 rays, etc. Figure 2 shows several 
families of rays (line drawings) for p-values ranging from one to 
eight. The value for the index of refraction is taken as 
n = 1.36 + Oi (i.e. an ethanol droplet) and the incident light is 
traveling from left to right in the figure. The line drawings of 
Fig. 2 show many different caustics. These drawings are two 
dimensional sections of the spherical droplet. The full ray 
picture is given by a figure of revolution about the Z axis for 
each drawing. Thus, the points where the p = 1 caustic intersects 
the droplet surface (on the shadow side) is a caustic ring when the 
figure is rotated. A photograph of the p = 1 caustic ring on an 
18.4 cm diameter thin-wall spherical glass globe filled with water 
is shown as Fig. 3a. The p = 2 caustic in Fig. 2 is a cusp of 
revolution whose focal point is on the Z axis at r = 0.61a when 
n = 1.36. A photograph of the p = 2 cusp of revolution for the 
thin-wall water-filled glass globe is shown as Fig. 3b. For each 
higher value of p, the caustic is also a cusp of revolution that 
winds further and further around the sphere, intersecting itself on 
the axis a number of times. These caustics will be referred to as 
the cusp caustics of the droplet. 

Employing the definition of a caustic as the envelope of the 
p-ray family, 2 we find in Appendix A that the coordinates of the p- 
ray cusp caustic are 

z = ( -sin (y) sin (0 r ) + k cos (y) cos (0 r ) ) a 

( 1 ) 

p = (-cos (y) sin (0 r ) - k sin (y) cos (0 r ) ) a 

where from Fig. 1, p is the distance perpendicular to the Z axis, 
6 r is the angle of refraction of the ray whose angle of incidence 
is 0 ; , 
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n sin (0 r ) = sin (0 i ) 


Y = Q d +pn - (2p-l) 0 r (3) 
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The cusp points of the caustics of Fig. 2 have the coordinates 


(-l) p a 
2p-l-n 


(5) 


Pc = 0 

and near a cusp point, the cusp takes the usual Pearcey form 3 


a 

where 

w = ^~ 1 ^ P ( 9 [ (2p-l -n) 3 - (2p-l-n) + n(n 2 -l) ] 

2 1 (2p-l-n) 4 

Since the caustics of Fig. 2 and eq. (1) correspond to p-rays 
transmitted into and reflecting within a single cross section of 
the sphere, the equations of the analogous ray theory caustics for 
a circular cylinder illuminated by normally incident rays are 
identical . 

In addition to the cusp caustic, a second type of caustic that 

is not readily apparent in Fig. 2, also exists within the droplet. 

It is the axial caustic. This caustic corresponds to focusing on 
the Z axis of light rays that enter the sphere along a ring of 
constant 6,. Thus for any p-family of rays, an axial caustic will 

be present if the p-rays cross the Z axis. For the droplet 

depicted in Fig. 2, it can be seen that p-values of two and larger 
have the required internal axis crossing which forms the axial 
caustic. Also note in Fig 2 that the axial caustic never extends 
across the full length of the droplet's axis since rays never cross 
the Z axis near the center of the droplet. 

The coordinates of the p-ray axial caustic within the sphere 

sin (0 ) 
z = -a — 

sin (y) (8) 

p = 0, 

are calculated in Appendix A. A photograph of the internal p = 2 
axial caustic for the thin-wall water-filled glass globe is shown 
as Fig. 3c. For = 0, the p-ray axial caustic begins at the cusp 
point of eq. (5) . For increasing 6, it proceeds along the axis 
inside the cusp of revolution toward the surface. For large p, it 
continues back and forth along the droplet axis with turning points 
at the locations where the cusp caustic intersects itself. This 
behavior is shown in Fig. 4 for the p = 5 axial caustic. In the 
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figure, the incident rays are labeled with the letters A-P, and the 
crossing point of each ray on the axis is labeled with the same 
letter. The exiting rays are denoted by the prime and are shown so 
that the direction of travel of the ray can be seen. Note in the 
figure that there are regions where the axis crossings by the p- 
rays are more dense: points A to D and K to N. It is within the 
second of these regions that the turning point (L) of the axial 
caustic occurs. Also note that some of the p = 5 rays do not cross 
the axis . These are rays that enter the droplet at angles greater 
than the entering angle of ray H and less than the entering angle 
of ray I; and also for rays that enter at angles greater than the 
entering angle of ray P. 

One important difference between the axial caustics and the 
cusp caustics is the degree of focusing that each represent. For 
the cusp caustics, incident rays that enter the droplet over a 
small element the droplet's surface area contribute to the caustic 
over a small element of the caustic's surface area. For the axial 
caustics, incident rays entering the droplet over an element of the 
droplet's surface area (in this case it is in the shape of an 
annular ring) contribute to the axial caustic over an element of 
the axial caustic's one-dimensional length. This mapping 
(surface -*• line) for the axial caustic corresponds to a higher 
degree of focusing than does the mapping (surface -*• surface) for 
the cusp caustics of revolution. As a result, the axial caustics 
are expected to be the largest contributor to the intensity on the 
droplet axis. The turning points of the axial caustics represent 
an even higher degree of focusing since at these points, incident 
rays entering the droplet over an annular ring are focused at the 
turning point rather than along a line. The method of Ref. 4 is 
employed in Appendix B to calculate the density of rays on the 
axial caustics. We find that the ray theory intensity of the p-ray 
axial caustics is given by 

dF p _ 2na 2 nKr(0 i ) i? p_1 (0j) sin (0 i ) cos (0 r ) sin 2 (y) C 9 ) 
dz cos (y) sin (0 r ) + k sin (y) cos (0 r ) 

where k is given in eq. (4) , and for unpolarized light T(0 ; ) and 
R(0j) are given by the averages of eqs. (18-38,39) of Ref. 1 
respectively . 

Wave Theory of the Caustics within a Liquid Droplet 

For a plane wave of wavelength X incident on a spherical 
droplet, the components of the interior electric and magnetic 
fields in Mie theory are given in Refs. 5-6,51-52. They are 
written as an infinite series of partial wave contributions which 
contain the Mie interior partial wave amplitudes 7 c t and d f . Each 
of the partial wave amplitudes in turn may be written as another 
infinite series corresponding to interactions with the droplet 
surface of the radially propagating spherical multipole partial 
waves. This series expansion of c, and d, is analogous to the 
situation for a flat thin film where the transmitted and reflected 
electric fields are decomposed into a series of terms corresponding 
to multiple internal reflections. 8 This series expansion for c t and 
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d f is known as the Debye series, and it provides the link between 
wave theory and the ray theory results for the p-ray families of 
Fig. 2. 910 The Debye series has been used for far-field scattering 
to isolate one component of the far-field intensity from 
interference effects produced between it and all the other 
components, and to suggest physical mechanisms for various 
scattering phenomena. 11 ' 20 It also provides a theoretical 
justification for the physical optics model for light scattering 
from bubbles. 21 ' 24 

Using the notation of Ref. 11, the Debye series expansions of 
the Mie theory partial wave amplitudes a f and b, for the scattered 
fields are 




1-Ri 
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( 10 ) 


and the Debye series expansions of the Mie theory partial wave 
amplitudes c e and d f for the interior fields are 
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1-R 
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In these expansions, region 1 is the interior of the sphere, 
region 2 is outside the sphere, T, 21 is the amplitude transmission 
coefficient for the i spherical multipole partial wave to be 
transmitted from region 2 into region 1, R, 11 is the amplitude 
transmission coefficient for the £ partial wave to be internally 
reflected from region 1 back into the same region, 


Tf = 
and 
R} 1 = 


-2i/.x 2 

(C(+tf) + i(tf-tf) 


- ( t| - ) - i ( t t + 1, ) 

.1 . x -2 \ 


(t i + tf) +i(tj 3 -t { 4 ) 

For a ( and c lr the t/ are given by 


t } = xy 


= xy 


t ( = xy 


J t (x) Jt-i(y) ~ nj t . x {x) j t (y) + j t ( x ) j, (y) 

n t (x) n t _ x (y) - hji m (x) n t (y) + n t (x) n t (y) 

n„(x) j,_i (y) - nn t _ x (x) j t (y) + n t (x) j t (y) 


j ( (x) n t _ x (y) - nj t . x (x) n t (y) + — n2 1) j t (x) n t (y) 


t, = xy 

and for b f and d f , the t/ are given by 


( 12 ) 


(13) 


( 14 ) 
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tf = xy[nn ( (x) n,^ 
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cf = xy [nj t (x) n t _ ± 
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„ - 2lza 

X ' 

y = nx, 


(y) - Ji-i (*) j, (y) ] 
(y) - n,.! (x) n,(y)] 
(y) - ^(x) j«(y)] 
(y) ~ jt-x(x) n t (y)] 


( 15 ) 


(16) 

(17) 


and where j f and n f are spherical Bessel functions and spherical 
Neumann functions respectively. (These are also known as spherical 
Bessel functions of the first and second kind.) In addition, the 
amplitude transmission and reflection coefficients satisfy 


l^ 22 l 

= i* ll i. 

(18) 

rpl2 _ 

if 1 

/ 

(19) 

n 


and 




<tt21 |2 

^-L + li ?, 11 ! 2 = i. (20) 


The physical interpretation of the various terms in eq. (10) are 
given in Ref. 11. For the interior of the droplet, the p term in 
eq. (11) corresponds to spherical multipole partial waves that have 
entered the droplet and have made p - 1 internal reflections. The 
Mie theory equivalent of the p-ray family of the previous section 
is the evaluation of the internal fields employing the p term in 
eq. (11) for c e and d f . 

Over the years, one concern in numerical Mie theory 
computations has been the determination of the partial wave £ max at 
which the Mie series is truncated. For far-field scattering, 

Dave's criterion 25 is that the Mie scattering amplitudes a f and 
satisfy |a fmax | 2 + |b fmax | 2 < 10' 14 . Occasional numerical difficulties 
with this criterion and the need for faster Mie calculation 
algorithms led Wiscombe to improve the criterion 26 to 
£ max = x + 4.05 x 1/3 + 2, which insures that a fmax and b fmax are each of 
order 10' 8 . For internal source function calculations, previous 
authors have used either the first 27 ' 28 or second 29 ' 30 criterion. 
However, each of these criteria leads to less precision for 
interior field calculations than for far-field scattering 
calculations. This can be seen from the Debye series as follows. 
For high partial waves (i.e., i > x) we have T f 21 -*• 0 and R, 11 -*• 1. 
Assume that for such a partial wave 

if 1 ~ e << l. (21) 


6 


Then substituting eq. (21) into eq. (20) and Taylor series 
expanding, we get 

i ?, 11 ~ 1 - — . ( 22 ) 

' 2n 

From eq. (11) the interior Mie amplitudes are of order 

C ( ~ d t ~ 6 , (23) 

while from eq. (10) the scattering Mie amplitudes are of order 

a, ~ b t ~ t 2 . (24) 

As a result, the Mie interior amplitudes which are calculated with 
c f and d f converge slower than the Mie scattering amplitudes which 
are calculated with a, and b f . To obtain the same precision for 
interior field computations that one demands for far-field 
scattering computations, more terms than indicated by Dave's or 
Wiscombe's criteria must be employed in the Mie series. A related 
situation occurs when searching for morphological scattering 
resonances in the far-field scattering amplitudes. 22,31 Using the 
formalism of Ref. 32, we find that the series cutoff 


<Lx = * + 7* 3 + 2 

ensures that c fmax and d fmax are of the order 10' 8 for 35 < x < 10,000 
and for Im(n) < 0.1. 

The cutoff of eq. (25) was used with the internal field 
formulas of Ref. 5-6 to determine the source function of a sphere 


5(f) = EUfyEjr) 

El 


(26) 


where E 0 is the incident field strength, for size parameter 
x = 100 . 0 (27) 

and real refractive index 


n = 1.36 (28) 

illuminated by unpolarized plane waves. The spherical Bessel 
functions that occur in the radial dependence of the electric field 
and in the evaluation of c t and d f were computed by downward 
recursion. The spherical Neumann functions were computed by upward 
recursion . 

The full Mie interior source function is shown in Fig. 5 in 
three different formats. Figure 5a is a 100 x 100 point carpet 
plot and shows the strong interior field enhancements on the 
sphere's axis at r « ± 0.8a. These enhancements in the source 
function are compatible with the locations of the enhancements 
found in previous Mie theory calculations 17 ' 33 ' 39 and will be addressed 
in the next section. Figure 5b is the same source function 
calculated on a 500 x 500 point grid and logarithmically mapped 
into 64 shades of gray: black indicating the most intense fields 
and white the least intense. This gray level plot emphasizes the 
interference structure of the interior source function. Figure 5c 
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shows the source function on the sphere axis. The enhancements at 
r ~ ± 0.8a appear to be a number of closely spaced individual 
enhancements . 

The contribution to the source function provided by the 
individual p = 1 to 8 terms of the Debye series of eq. (11) were 
also calculated for unpolarized incident light and are shown as the 
gray level plots in Fig. 2 along with the ray theory predictions. 
A common rule of thumb has been that ray theory becomes a good 
approximation to the far-field primary rainbow caustic for 
x > 2000. 40 Furthermore, when rainbows of higher order are 
observed, larger values of x are required before these rainbows are 
well formed. 25 For the interior caustic structure of Fig. 2, ray 
theory is already a qualitatively good approximation to wave theory 
at x = 100 for all values of p up to at least 8. 

Fig. 2 shows a number of other interesting features. As with 
all other wave theory optical caustics, 41 the infinite intensity of 
ray theory caustics is rendered finite and is surrounded by 
diffraction and interference structure. For p = 2, the interior 
portion of the first order rainbow is visible in both the ray 
picture and the Debye series gray level plot. It begins somewhat 
above and below the p = 2 cusp point and extends to the sphere 
surface, where it continues into the far-field to become the 
familiar first order rainbow. This behavior is shown more visibly 
in the p = 2 carpet plot of Fig. 6. Also surprisingly visible in 
Fig. 6 are the first few supernumerary rainbows inside the droplet. 
One other feature to note in Fig. 2 is how the interference 
structure in the gray level plots corresponds to regions of 
intersecting rays in the line plots. For example, the gray level 
plot for p = 3 shows much interference structure on the illuminated 
side of droplet (from r « -1.0a to r « -0.27a) and a smooth 
decrease in intensity on the shadow side (from r ~ -0.27a to 
r « +1.0a) . The corresponding line plot for p = 3 shows many rays 
intersecting one another in the region on the illuminated side and 
the rays diverging in the region on the shadow side. 

In order to assess the quantitative agreement between the 
caustics of ray theory and the individual terms of the Debye series 
expansion of the interior electric field, the source function was 
calculated along the Z axis for the p = 1 to 10 terms of the Debye 
series. The results are shown in Fig. 7 as the solid curves. It 
is remarkable that for large p, the source function for an 
individual term in the Debye series, ranges over as many as sixteen 
orders of magnitude. This large variation is easily understood. 
From the previous section, the intensity of a p-ray is proportional 
to T ( 0|) R p_1 ( 0j) . Consider the refractive index of eq. (28) and 
p = 10. The intensity at the center of the sphere is due to rays 
incident near the sphere axis with 6 { ~ 0, T ~ 0.977 , R ~ 0.023, and 
TR 9 « 10 15 . But the intensity near the surface of the sphere is due 
to rays near grazing incidence with ~ n/2, and perhaps T « 0.01 
and R « 0.99 giving TR 9 ~ 10' 2 , thirteen orders of magnitude larger 
than at the center. In other words, rays near grazing incidence 
have a low probability of being transmitted into the droplet. But 
once inside, the probability of their surviving many internal 
reflections is quite large. On the other hand, rays near normal 
incidence are very likely to be transmitted into the droplet. But 


8 


once inside, they are again very likely to be transmitted back out. 
Since the cusp points of the p-ray caustics are due to rays near 
normal incidence, this argument explains why they are the locations 
on the caustics of nearly the weakest intensity, rather than the 
strongest. The illusion in the ray theory pictures of Fig. 2 that 
the cusp point is the most intense location on the caustic results 
from incorrectly drawing all the ray lines with equal darkness or 
thickness. Since T (0 t ) R pl ( 0j) continues to grow in going from normal 
incidence to grazing incidence, the axial caustics become more 
intense as they progress from the center of the droplet to the 
perimeter . 

In the previous section, we argued that on the Z axis, the 
electric field should be dominated by the contribution provided by 
the axial caustics. Figure 7 shows the density of rays on the 
axial caustics (data displayed as points) given by eq. (9) 
multiplied by a single scaling factor for all p . The locations on 
the curves that are multivalued represent regions where the axial 
caustic has turned back on itself. As is seen in this figure, at 
x = 100, the ray theory model of caustics within a liquid droplet 
is a quantitatively good approximation to the caustic structure 
predicted by wave theory. At the turning points of the axial 
caustics, the ray theory intensity becomes infinite while the Debye 
series intensity undergoes a rapid increase by a number of orders 
of magnitude and is accompanied on one side by an interference 
structure, a behavior typical of diffraction caustics. 41 For 
larger x the interference structure becomes finer and the caustics 
become quite well defined. 

Comparison with SRS and Fluorescence Emission 

The transition rates for stimulated Raman scattering (SRS) and 
fluorescence emission at a given location within a sample are 
proportional to the square of the stimulating electric field at 
that location, i.e. the source function. 42 Exploiting this fact, 
these inelastic light scattering processes have been used to probe 
the distribution of electric field strength within spherical liquid 
droplets illuminated by laser light. 33 ' 344344 Three regions of 
enhanced source function within the droplet have been observed. 
The first region is on the droplet's shadow side at r « +0.8a. The 
second region is on the droplet's illuminated side at r « -0.8a. 
The third region is the entire circumference of the droplet. This 
last region of strong emission, in contrast to the first two, is 
observed only for certain diameter droplets and corresponds to the 
situation where the diameter of the droplet, its refractive index, 
and the wavelength of the incident light satisfy the conditions for 
the occurrence of a morphological scattering resonance. 4547 

The analysis of the previous two sections and Fig. 7 can be 
used to explain the first region of enhanced source function on the 
droplets shadow side. The regions of most intense focusing on the 
shadow side of the droplet axis in Fig. 7 occur at focal point of 
the p = 2 rays, and also at the maxima of the Debye series 
interference patterns associated with the outermost turning points 
of the axial caustics at r a +0.80a for p = 6 and r » +0.78a for 
p = 10. Other turning points in this .region are the outermost 
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turning points of the p = 5 and 9 axial caustics. But these are 
lower than those for the p = 6 and 10 by about two to five orders 
of magnitude. As a result we attribute the r ~ +0.8a source of SRS 
and fluorescence emission to the p = 2 focal point amplified by the 
near-identical turning points of the p = 6, 10, 14, ... axial 
caustics . 

In Refs. 33-34, 43 the shadow side source of SRS and 
fluorescence emission seems to be on the droplet axis. However in 
Ref. 44 it appears instead on the surface of the droplet at the 
position of the p = 1 caustic ring . Kerker and his co-workers 
have pointed out that the electric and magnetic field boundary 
conditions at the surface of the droplet must be satisfied not only 
by the radiation at the incident wavelength, but also by the 
radiation at the inelastically scattered wavelength. 42,48 ^ 19 This 
causes the inelastically scattered radiation to be preferentially 
emitted in certain directions, and that these directions vary for 
stimulated molecules at different locations within the droplet. 
Since the observations of Ref. 33-34,43 and Ref. 44 were made at 
different scattering angles, this directionality of the 
inelastically scattered radiation may provide the explanation for 
the difference between the observations. This point, however, 
requires further study. 

The second region of SRS and fluorescence emission on the 
illuminated side of the droplet (at r ~-0.8a) is said to be caused 
by the focusing on the axis of rays that have internally reflected 
twice within the droplet 35 (i.e. p = 3 rays). This interpretation 
can be seen in Figs. 2, 7 to be incorrect since for n = 1.36, eq. 
(5) predicts that the p = 3 cusp point should occur at z = -0.27a. 
This corresponds to the very small enhancement on the full Mie plot 
in Fig. 5c at z ~ -0.27a. The locations of the outermost turning 
points of the axial caustics at r ~ -0.80a for p = 4 and at 
r ss -0.78a for p = 8 closely correspond to the maximum source 
function enhancement in Fig. 5c. As a result we attribute the 
r « -0.8a source to the nearly identical locations of the turning 
points of the p = 4, 8, 12 ... axial caustics. 

It has been noted that when observing the sources of SRS and 
fluorescence emission within liquid droplets, the curvature of the 
droplet surface causes the apparent source of the emission to be 
shifted from its actual location. 33,50 This shift in position is 
calculated in Appendix C using geometrical optics and is shown in 
Fig. 8 for n = 1.333. Also shown in Fig. 8 are the results of an 
experimental verification of the position shift. The 18.4 cm 
diameter thin-wall glass globe described earlier was fitted with a 
scale along its diameter. The globe was photographed from the 
side, then filled with water and photographed again. The shifts in 
the positions of the grid marks on the scale between the two 
photographs were measured. The agreement with the results of 
Appendix C is good. 

We applied the method of Appendix C to determine the actual 
locations of the sources of SRS and fluorescence emission in the 
photographic observations of Refs. 33-34, 43. The results are 
given in Table 1. The experimental results agree with the 
locations of the p = 2 cusp point and the outermost turning points 
of the p = 4, 6, 8, 10, ... axial caustics found in Fig. 7. 
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Conclusions 


In the past, the Debye series expansion of the far-field Mie 
scattering amplitudes has provided a valuable method by which a 
scattering mechanism can be computationally isolated from all the 
other effects that interfere with it. In this paper we show that 
the Debye series expansion of the interior field Mie amplitudes has 
a physical interpretation of equal richness. The individual terms 
of the Debye series correspond to families of rays inside the 
droplet. The caustics of ray theory are reproduced by the Debye 
series calculations and decorated with diffraction structure 
typical of wave theory treatments. Further, the Debye series 
analysis provides the physical explanation of the source function 
hot spots observed using inelastic light scattering. 
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Appendix A 


(Al) 


Consider the p-ray which enters the liquid droplet with the 
angle of incidence 6 { at the coordinates 

z ± = - a cos (0 i ) 

Pi = +<a sin (0 i ) . 

The location where this ray intersects the sphere surface at the 
p - 1 internal reflection is 

z p -x = ~a c os(0 i + 2(p-l)4>) 

Pp_i = + a sin (0 i + 2(p-l)<j>) 

and the location where it intersects the sphere surface at the p 
internal reflection is 


(A2) 


z p = -a cos (0i + 2 p<J>) 

P p = +<a sin(0 i + 2 p4>) , 
where 


(A3) 


4 > = § - e r (A4) 

and 6 r is the angle of refraction given in eq. (2). The equation 
of this p-ray is 


P = Pp-i - [tan (y) ] {z-z p _x), (A5) 

where y is given in eq. (3) . Consider the intersection of this ray 
with the p-ray whose angle of incidence is 0 , + e. The angle of 
refraction of this second p-ray is 6 t + S where 


e cos (©j) 
n cos (0 r ) 


(A6) 


for small e. In the limit as e -*• 0, the intersection point of the 
first and second p-ray is the location on the caustic given by 
eq. (1). 

For the axial caustic, consider the location where the p-ray 
crosses the Z axis. Eq. (A5) with p = 0 gives 


z = 


z p - i + 


j p - i 


tan (y) 
Substituting eq. (A2) 


into (A7) gives eq. (8) . 


(A7) 
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Appendix B 


Let dF^ be the flux incident on the annular ring of width d0; 
centered at the angle of incidence 9 t and let dz be the progress 
along the p-ray axial caustic that results from increasing the 
angle of incidence from 9 { - d9j2 to 9, + d9 t /2. Then 

dF inc = 2 n a 2 sin(0 i ) cos (0 i ) d0 i . (Bl) 

The fraction of this incident flux dF p that survives p-1 internal 
reflections is 


dF p = dF lnc T(6 1 ) RF-iiQJ . (B2 ) 

Then following the procedure of Ref. 4, 


dF^ m dQ, 
dz dz 

dd, 


(B3) 


The expression for z is given in eq. (8) . Differentiating it with 
respect to 9, and inserting the result into eq. (B3) gives eq. (9) . 
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Appendix C 


Consider a point on the axis of the sphere at the location 

z s = wa (Cl) 

as in Fig. 9a. Assume that this point source emits light rays in 
all directions. Consider the light ray that leaves the source at 
an angle a with respect to the p axis. Then the z coordinate of 
the intersection of this light ray with the surface of the sphere 
is 


where 

u = w + [w 2 + 1 -2w cos (P) ] 1/4 cos (P + 0 i ) 
and 

sin<e,) = 51sin(P) 

[ w 2 + 1 - 2 w cos ( P ) ] ' /s 

The angle that the exiting ray makes with the Z axis is 

l = P + 9r 

where 


sin (0 r ) 


n w sin (P) 

[ w 2 + 1-2 wcos (P) ] ' h 


(C2 ) 
(C3 ) 

(C4 ) 

(C5) 

(C6) 


As is shown in Fig. 9b, many such light rays leave the source S, 
pass through a lens. These rays are then imaged at the point S on 
a viewing screen parallel to the Z axis as shown in Fig 9c. These 
rays, such as rays A, B, and C in Fig. 9b appear to originate from 
the virtual source S' within the droplet. Consider ray B which 
passes through the center of the lens. If d is the distance from 
the droplet to the lens, then from Fig. 9c the angle that this ray 
makes with the Z axis is approximately 

5 - -§ + arctan(^). (C7) 

In the limit d >> a > z out , combining eqs. (C5, C7) gives 
cos(P) = sin(0 r ) (C8) 

or 

w 2 + 1 - 2 w cos ( P ) = n 2 w 2 tan 2 ( P ) . ( C9 ) 

Finally, given the actual position of the source w, one numerically 
solves eq. (C9) to obtain B. This is then inserted into eq. (C4) 
to find 6 ;, and the result is substituted into eq. ( 03 ) to find u. 
If a << d, this is also the apparent position of the virtual 
source S ' . 
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Table 1 


Ref. 

Liquid 

Droplet 

Radius 

Shadow Side 
Enhancement 

Illuminated Side 
Enhancement 

33 

ethanol 
n=l . 36 

35 /jm 

0.66a < r < 0.79a 

-0.78a < r < -0.70a 

34 

methanol 
n=l . 33 

25 fim 

0.67a < r < 1.00a 

-1.00a < r < -0.63a 

43 

ethanol (4) 
n=l . 36 

not 

given 

0.42a < r < 1.00a 

-1.00a < r < -0.26a 


+ The experimentally observed locations of the sources of SRS 

and fluorescence emission in liquid droplets. 

(a) P. Chylek, private communication. 



Figure 1 A geometrical light ray with an angle of incidence 0 i 
entering and internally reflecting within a liquid 
droplet of radius a. 
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Figure 2a - Gray level plot of the p = 1 term of the Debye Series 
interior source function within the droplet for x = 100 
and n = 1.36, for unpolarized light. 
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Figure 2b - Line Drawing of the geometrical p = 1 ray family within 
the droplet for n = 1.36. 



Figure 2c - Gray level plot of the p = 2 term of the Debye Series 
interior source function within the droplet for x = 100 
and n = 1.36, for unpolarized light. 
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Figure 2d - Line Drawing of the geometrical p = 2 ray family within 
the droplet for n = 1.36. 
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Figure 2e - Gray level plot of the p = 3 term of the Debye Series 
interior source function within the droplet for x = 100 
and n = 1.36, for unpolarized light. 
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1 



Figure 2f - Line Drawing of the geometrical p = 3 ray family within 
the droplet for n = 1.36. 
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Figure 2g - Gray level plot of the p = 4 term of the Debye Series 
interior source function within the droplet for x = 100 
and n = 1.36, for unpolarized light. 
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Figure 2h - Line Drawing of the geometrical p = 4 ray family within 
the droplet for n = 1.36. 
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Figure 2i - Gray level plot of the p = 5 term of the Debye Series 
interior source function within the droplet for x = 100 
and n = 1.36, for unpolarized light. 
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Figure 2j - Line Drawing of the geometrical p = 5 ray family within 
the droplet for n = 1.36. 


29 





Figure 2k - Gray level plot of the p = 6 term of the Debye Series 
interior source function within the droplet for x = 100 
and n = 1.36, for unpolarized light. 


30 



ipg 

SOcrtf 

■ mm* — 

m 

1 

fill 

■ 

Mm 


•-V y 

U !»< 

f.tM'A 

Pm 



9 


JVi 

K! 


Figure 2L - Line Drawing of the geometrical p = 6 ray family within 
the droplet for n = 1.36. 
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Figure 2m - Gray level plot of the p = 7 term of the Debye Series 
interior source function within the droplet for x = 100 
and n = 1.36, for unpolarized light. 



Figure 2n - Line Drawing of the geometrical p = 7 ray family within 
the droplet for n = 1.36. 
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Figure 2o - Gray level plot of the p = 8 term of the Debye Series 
interior source function within the droplet for x = 100 
and n = 1.36, for unpolarized light. 
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Figure 2p - Line Drawing of the geometrical p = 8 ray family within 
the droplet for n = 1.36. 
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Figure 3 (a) The p = 1 caustic ring, (b) the p = 2 cusp caustic 

of revolution, and (c) the p = 2 axial caustic for an 
18.4 cm thin-wall water-filled spherical glass globe. 
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Figure 4 Formation of the p = 5 axial caustic for various incident 
ray angles 0^ Turning points occur at the locations 
where the p = 5 cusp caustic intersects itself on the Z 
axis. The arrows below the figure show the progression 
of the axial caustic for the rays A-H and I-P. 
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Normalized Distance (r^d) 


Figure 5 The interior source function in Mie theory for a sphere 
with x = 100, n = 1.36, and for unpolarized incident 
light shown (a) as a carpet plot, (b) as a gray-level 
plot, and (c) along the Z axis of the droplet. 
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Figure 6 The p = 2 Debye series contribution to the interior 

source function showing the p = 2 cusp point and the 
interior portion of the first order rainbow and 
supernumerary rainbows. 
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10 










Figure 7 - 

The contribution of the p-term of 
the Debye series to the interior 
source function on the droplet 
axis for p = 2 to 10 (curve) and 
the density of geometrical rays 
on the axial caustic from 
eg. ( 9 ) (points) . 
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Apparent Position (r/a) 



Figure 8 The apparent position of a light source on the droplet 
axis as a function of the actual position for the 
geometric optics model of Appendix C and n = 1.33. The 
solid circle data points are the measured positions of 
grid marks on a scale fitted along the diameter of an 
18.4 cm thin-wall glass globe when empty (actual 
positions) and filled with water (apparent positions) . 
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Figure 9 





(a) A light ray leaving a source at the coordinate z s on 
the axis of a spherical droplet. (b) Three such rays 
leaving the source S and imaged by a lens. The virtual 
source of the rays A, B, and C is S'. (c) The imaging of 
the spherical droplet by a lens. _ The image of the 
topmost portion of the droplet T is T. The image of the 
lowermost portion of the droplet L is L. The image of 
the source S is S. 
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